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We investigate two-pion Bose-Einstein correlations of quark-glnon plasma droplet sources in single¬ 
event measurements. We find that the distribution of the fluctuation between correlation functions 
of the single- and mixed-events provide useful signals to detect the granular structure of the source. 
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I. INTRODUCTION 

Recently, there has been much progress in the exper¬ 
imental search for the quark-gluon plasma [IQ While 
there are many indications suggesting the presence of a 
very dense matter produced in high-energy heavy-ion col¬ 
lisions QQ, a very important question is whether the 
produced dense matter is a quark-gluon plasma. If it 
is a quark-gluon plasma, it will undergo a phase transi¬ 
tion from the quark-gluon plasma phase to the hadronic 
phase. It is desirable to search for the signature for the 
phase transition of the quark-gluon plasma. 

The signature for the phase transition however depends 
sensitively on the order of the transition. Previously, it 
was suggested by Witten and many other workers that a 
granular structure of droplets occurs in a first-order QCD 
phase transition, and the observation of the granular 
structure can be used as a signature for a first-order OCD 
phase transition Q,QQIEQ,0,0,0,0)1313 QlS- 

In a recent analysis, a granular emitting source of 
droplets was put forth to explain the HBT interferometry 
puzzle for nucleus-nucleus collisions at RHIC ^5. The 
suggestion was based on the observation that in the hy- 
drodynamical model 013, the particle emission time 
scales with the radius of the droplet. Particles will be 
emitted earlier if the radius of the droplet is smaller, as 
in a source of many droplets. An earlier emission time 
will lead to a smaller extracted HBT radius i?out- As 
a result, the value of Rout can be close to i?side for a 
granular quark-gluon plasma source |l 4| . 

Previously, Pratt et al. studied the HBT interferome- 
t^ of granular droplets by averaging over many events 
0 . Methods to detect a granular structure by the single¬ 
event intensity interferometry were recently proposed 
0 It was found that the single-event correlation func¬ 
tion from a chaotic source of granular droplets exhibits 
large fluctuations, with maxima and minima at relative 
momenta which depend on the relative coordinates of the 
droplet centers. The presence of this type of maxima and 
minima of a single-event correlation function at many rel¬ 
ative momenta is a signature for a granular structure and 
a first-order QCD phase transition 0 . 

The difficulty of using the single-event two-pion inter¬ 


ferometry at RHIC arises because of the small number 
of pion pairs with small relative momenta. In a typical 
single event of a nearly head-on collision at very high en¬ 
ergies at RHIC, the number of identical pions is of the 
order of a few thousand. The number of observed identi¬ 
cal pions tItt is only a small fraction of this number. For 
example, the number of identical pions detected in the 
STAR Collaboration in the most central Au-Au collisions 
at RHIC is of the order of a few hundred 0 . Although 
the number of pairs of identical pions in the event varies 
as — l)/2, only a small fraction of these 

pairs have relative momentum small enough to be useful 
in a HBT analysis. The number of pion pairs in each 
relative momentum bin may be so small that there can 
be large associated statistical errors. 

Instead of trying to obtain the detailed granular struc¬ 
ture of the emitting source in each event at present, it 
will be useful in the initial stage to have a more modest 
goal. It is desirable to see whether the correlation func¬ 
tions indicate possible signals for a granular structure. 
More refined study of the granular structure can follow 
after the initial stage becomes successful. 

Accordingly, we shall try to outline a method to ex¬ 
tract the signals for the granular structure. Our idea is 
to calculate the fluctuations of the single-event correla¬ 
tion function relative to its corresponding mixed-event 
correlation function. The difference constitutes the “sig¬ 
nal” for the event in question. The distribution of these 
fluctuations (in units of their statistical errors), collected 
for a large number of single-events to enhance statistics, 
would have a wider distribution for a granular struc¬ 
ture, compared to those from emitting source without 
the granular structure. The distribution of the correla¬ 
tion function fluctuations provide a useful tool to detect 
the granular structure of the source. In particular, the 
root-mean-square fluctuation, in the case of a granular 
droplet structure, increases when the number of droplets 
decreases. If the phase transition is accompanied by only 
a few number of droplets, the signal may be large enough 
to make it detectable. In what follows, we would like to 
analyze whether this method may be feasible. 
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II. SINGLE-EVENT AND MIXED-EVENT 
TWO-PION CORRELATION FUNCTIONS 


The two-particle Bose-Einstein correlation function for 
the detection of identical pions with laboratory momenta 
ki and k 2 is defined as C{ki,k 2 ) = P{ki, k 2 )/P{ki)P{k 2 ) 
where P{ki,k 2 ) is the two-particle momentum distribu¬ 
tion and P{ki) is the single-particle momentum distri¬ 
bution with momentum ki. For a chaotic pion-emitting 
source, P{ki) {i = 1,2) is 

P{h)=Y,A\h,X,), (1) 

Xi 


where A{ki,Xi) is the magnitude of the amplitude for 
emitting a pion with 4-momentum ki = (Ei^ki) at Xi. 
The two-particle distribution function P(ki,k 2 ) can be 
expressed as 


P(fcl,fc2)= ^ |$(fci,fc2;Xi,A2) 

Xl,X2 


2 

( 2 ) 


where ^{ki, k 2 ', Xi, X 2 ) is the two-pion wave function. 
Neglecting the absorption of the emitted pions by other 
droplets, ^{ki, k 2 ] Xi, X 2 ) is simply 


A:2; ^ 1 , -^ 2 ) 

= ^ \Aiki,Xi)A{k2, 

V 2 L 

-f A{ki,X2)A{k2,Xi)P’^^-^^+^'^^-^^ 


( 3 ) 


The correlation function C{ki, ^ 2 ) is in general a func¬ 
tion of the 4-dimensional momentum ki — /c 2 and ki + k 2 . 
Previous results of the single-event correlation function 
of granular droplets indicate large fluctuations as a func¬ 
tion of the relative 4-momentum ki — k 2 , having max¬ 
ima and minima at locations which depend on the rel¬ 
ative coordinates of the droplet centers El To map 
out the details of such a multi-dimensional correlation 
function, it is necessary to have a large number of pion 
pairs in a single event which may be beyond the capabil¬ 
ities of present detectors and accelerators. With limited 
statistics as would likely be the case, we can first study 
the simplifying case by concentrating on a small number 
of degree of freedom and integrating out other degrees 
of freedom so that the statistical errors in the correla¬ 
tion function can be smaller. For this purpose, we shall 
study the correlation function as a function of the vari¬ 
able q = |ki — k 2 | for which the other degrees of freedom 
have been integrated out. 

In our numerical work for a granular source with 
droplets, we obtain the single-event and mixed-event two- 
pion correlation functions for the granular source with 
the following steps: 

Step 1: Generate the space-time coordinates Rj {j = 
1 , 2 ,..., Nd) of the droplet centers according to a distri¬ 
bution of the droplet centers. 


Step 2: Select the two emitting-pions from the droplets 
randomly, and get the space-time coordinates Xi and 
X 2 of the pions in the laboratory frame according to the 
density distribution of the droplets, the coordinates of 
the droplet centers, and the collective velocities at the 
emission points v{Xi) [i = 1,2). 

Step 3: Generate the momenta k'^ = (if), k)^) and /c) = 
(if), k)) of the two pions in the frame in which the source 
element is at rest according to the distribution A^{k[,X), 
taken to be the Bose-Einstein distribution characterized 
by the temperature Ty, and obtain their momenta ki and 
k 2 in the laboratory frame by Lorentz transformation if 
the source element is boosted in the laboratory frame. 

Step 4: Accumulate the event in the bin of the corre¬ 
sponding relative momentum variable q = |ki — k 2 | with 
the weight factor wn for the probability P{ki)P{k 2 ) for 
a pair of uncorrelated pions of relative momentum q, and 
with the weight factor W 12 for the probability P{ki,k 2 ) 
for a pair of correlated pions of relative momentum q, 

u;n = [if)/ifi][if)/if 2 ], (4) 


W,2 = 2 


Aki-Xi-\-ik2 •X2 


^E'JE^^E'JE2 h 
+ ^E[.2/Ei^E'.2i/E2 


( 5 ) 


where E[ - {i,j = 1, 2) is the energy of the i—th pion in the 
frame in which the source element at Xj is at rest, which 
is obtained from Ei by a reverse Lorentz transformation 
with the collective velocity v{Xj). 

Step 5: Repeat steps 2 through 4 for pairs of pions 
in a single event. We label the distributions obtained for 
P{ki)P(k 2 ) for a pair of uncorrelated pions by Uncors(( 3 ') 
and the distribution of P(A:i, fe) for a pair of correlated 
pions by Cors{q). 

Step 6: Repeat steps 1 through 5 for A"event number 
of different events, and obtain the mixed-event corre¬ 
lated and uncorrelated pion-pair distributions CoTm{q) 
and Uncorm(9), by summing Gors(g) and Uncors(( 7 ) of 
the Vevent different events. 

Step 7: Obtain the single-event correlation function 
Cs{q) by dividing Cois{q) by Uncorm(9) and the mixed- 
event correlation function Cm{q) by dividing CoTm{q) by 
Uncorm(g). 

We first investigate the two-pion correlation functions 
for Nd static granular droplet sources. The centers of the 
droplets are assumed to follow a Gaussian distribution 
with a standard deviation cr/j, and the density distribu¬ 
tion of each droplet is assumed to be given by a Gaussian 
distribution with a standard deviation ad. In our numer¬ 
ical examples, an and ad are taken as 5.0 and 1.5 fm, and 
the thermal emission temperature of the pions is taken 
to be O.GbTc = 0.65 x 160 = 104 MeV. Fig. 1 (a), (b), 
and (c) show the two-pion correlation functions for the 
granular sources with Nd =4, 8, and 16, respectively. In 
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FIG. 1: Two-pion correlation functions for a sample of differ¬ 
ent single events (dashed lines) and mixed events (solid lines) 
for static granular sources of Nd droplets. 



/ 

FIG. 2: The distributions dN/df of / for (a) = 5 x 10^ 

and (b) Nt„t = 10®. 


III. DISTRIBUTION OF THE FLUCTUATION 
OF SINGLE-EVENT TWO-PION CORRELATION 
FUNCTION 

Form the results in Fig. 1, we observe that the cor¬ 
relation function Cs{q) for individual single events fluc¬ 
tuates with respect to the mixed-event correlation func¬ 
tion Cm{q)- We can make the fluctuation quantitative 
and introduce the fluctuation as the difference between 
the single-event correlation function and the mixed-event 
correlation function. In order to take into account the er¬ 
ror of the measurement, we weigh the fluctuation by the 
inverse of the corresponding error and define the fluctu¬ 
ation quantitatively as 


\Cs{q)-C^{q)\ 

^\Cs{q)-Cm{q)\' 


( 6 ) 


where A\Cs{q) — Cm{q)\ is the error in the measurement 
of Cs{q) - Cm{q) given by 


each figure the dashed lines give the correlation function 
C(q) for a sample of different single events each of which 
has a correlated pion-pair distribution Cors(( 7 ) calculated 
with Ntttt = 10® pion pairs within q < 250 MeV, and the 
solid line is for the mixed-event obtained by averaging 
Nevent = 10® single events. One can see that there are 
fluctuations for the single-event correlation functions rel¬ 
ative to the mixed-event correlation function, and the 
fluctuations increase as Nj, decreases. 


A|a(g) - Cm{q)\ = ACs{q) + ACm{q) ^ ACs{q) 

We calculate the distribution of the fluctuation / for 
the sources with different numbers of droplets. In these 
calculations, we take the width of the relative momentum 
bin as 10 MeV and use the bins in the region 20 < q < 250 
MeV. Fig. 2 (a) and (b) show the distributions of / for 
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FIG. 3: The root-mean-square / as a function of Nd for static 
granular sources. The double-dot-dashed lines are the result 
for a single Gaussian source. 


Nd = 4, 8,16, and 32, obtained from 1000 single events 
each of which has the correlated pion-pair distribution 
Cors(g') calculated with = 5 x 10"* and 10® pion 
pairs within q < 250 MeV. The standard deviations an 
and Gd for the granular sources are 5.0 and 1.5 fm. The 
results for a non-granular single source with a Gaussian 
density distribution with 5.0 fm standard deviation are 
also shown as a reference. It can be seen that the distri¬ 
butions for the granular sources are wider than that for 
the single source. The width of the distribution for the 
granular source decreases with Nd- 

From the distribution of /, one can calculate the root- 
mean-square values of /. Fig. 3 (a) and (b) show the 
root-mean-square / as a function of Nd for sources with 
CTfl =5.0 and 7.0 fm obtained for for the cases of N^^^ = 
10® and ad =1.5 fm, and Nt^.,^ = 5x 10"^, and ad =2.0, 1.5, 
and 1.0 fm. The double-dot-dashed lines are the root- 
mean-square / for the distribution of the single Gaussian 
source in Fig. 2. It can be seen that /rms is sensitive to 
Nd and decreases as ad and a^ increase. With the values 
of /rms, one may distinguish the granular source with 
Nd=l& from the single Gaussian source with iV,rir = 5 x 
10^ pairs of identical pions, and one may even distinguish 
the granular sources up to Nd = 32 with iV,r 7 i- = 10® 
identical pions. 


IV. DISTRIBUTION OF THE FLUCTUATION 
FOR HYDRODYNAMICAL QGP DROPLET 
SOURCE 

We investigate next the distribution of / for a granu¬ 
lar sources of quark-gluon plasma droplets which evolve 
hydrodynamically. We assume that all of the droplets in 
the source have the same initial radius Vd = 1.5 fm and 
evolve hydrodynamically in the same way. We use rela¬ 
tivistic hydrodynamics with the equation of state of the 
entropy density 1^ to describe the evolution of the 
droplets [ij, [l^ [ 13 , and take the temperature width of 
the transition as AT = 0.05 and the initial conditions of 
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FIG. 4: Two-pion correlation functions for a sample of differ¬ 
ent single-events (dashed lines) and mixed events (solid lines) 
of the dynamical granular sources. 


the droplets as [Hill HI 




eo, r' < rd, 
0, r' > rd, 


v'{t' = 0,r') 


0, r' < rd, 

1, r' > rd, 


( 8 ) 


where, r', e', and v' are the radial coordinate, energy 
density, and velocity of a fluid element in the droplet- 
rest frame, eo = 1.875rcSc [3113 is the initial energy 
density of the droplets. The initial distribution of the 
droplet centers is taken to be a Gaussian distribution 
with the standard deviation an = 5.0 fm. For the case 
with an additional collective radial expansion, the droplet 
centers are assumed to have a constant radial velocity Vd 
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in the center-of-mass frame of the granular source |14| . 
In order to reduce the influence of source lifetime on our 
observations for the granular structure we use the 
“side” component of the relative momentum of the two 
pions Qs (perpendicular to the total momentum of the two 
pions) as the variable for the granular 

source of hydrodynamic evolution QGP droplets. Fig. 
4 (a), (b), and (c) show the correlation function C{qs) 
(with gout < 20 MeV) for the dynamical granular sources 
with Vd = 0.5 and Nd =4, 8, and 16, respectively. In each 
figure the dashed lines give the correlation function C{qs) 
for a sample of different single events each of which has a 
correlated pion-pair distribution Cors(gs) calculated with 
= 10® pion pairs within qs < 250 MeV, and the 
solid line is for the mixed-event obtained by averaging 
Nevent = 10^ single events. The freeze-out temperature 
of the pions is taken to be O.GSTc = 0.65 x 160 = 104 
MeV. 


Fig. 5 (a) and (b) show the distributions of / for the 
dynamic granular sources with Vd = 0 and Vd = 0.5, and 
Nd = 4, TVd = 8, Nd = 16, and Nd = 32. The double- 
dot-dashed lines are for a dynamical single source with 
initial radius rd =5.0 fm as a reference. The number of 
pion pairs within qs < 250 MeV for one single event is 
Ntttt = 5 X 10^ and the number of events is Vevent = 10^ 
for all the distributions. Because of the radial expan¬ 
sions of the droplets, the two pions emitted from the 
single-droplet source are boosted along the radial direc¬ 
tion of the source and the two pions emitted from the 
granular source are boosted along random direction when 
Vd = 0.5. It leads to a smaller width of the distribution 
for the single-droplet source than those for the granular 
sources. When Vd 0 the widths of the distributions for 
the granular sources decrease because of the additional 
collective radial expansion Vd- 


Fig. 6 (a) and (b) show the root-mean-square val¬ 
ues of / calculated from the distributions of one thou¬ 
sand single-events for the dynamical granular sources 
with Vd =0 and 0.5, respectively. The symbols of A are 
for the case of the number Nt^-tt of the pion-pair within 
qs < 250 MeV for one single event is 10®. The sym¬ 
bols of • and * are for the cases of the N^^tt = 5 x 10"'^ 
and with and without the constraint qo < 20 MeV. The 
double-dot-dashed lines are the result for the distribu¬ 
tion of the single-droplet source in Fig. 5. One observes 
that /rms for the granular sources are larger than that 
for the single-droplet source, and the root-mean-square 
/ for the granular sources increase with Although 

the /rms results with the constraint qo < 20 MeV are 
larger than those without the constraint, our simulations 
indicate that this constraint is very restrictive and only 
about 3.5 percent of the pion pairs within qs < 250 MeV 
satisfy it. As the droplet number decreases, the fluctu¬ 
ation increases. The effect of the increase becomes less 
pronounced as the collective expansion velocity increases. 
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FIG. 5: The distributions of / for the dynamical quark-gluon 
plasma granular sources with Vd =0 and 0.5. 
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FIG. 6: The root-mean-square / as a function of Nd for 
dynamical granular sources. 


V. CONCLUSIONS AND DISCUSSIONS 


Recent experiments at RHIC provide ample evidence 
for a dense matter produced in high-energy heavy ion 
collisions. Is the produced dense matter the quark-gluon 
plasma? If so, what is the order of its phase transition? 
As a granular structure of droplets occurs in a first-order 
QCD phase transition, the observation of the granular 
structure can be used as a signature for a first-order OCD 
phase transition i i i i i i m HI in in niiii 

We would like to develop tools to use HBT interfer¬ 
ometry to examine the granular droplet structure of the 
dense matter if they are produced in high-energy heavy- 
ion collisions. We show previously that the single-event 
correlation function in HBT interferometry for granular 
droplets exhibits oscillations, depending on the relative 
coordinates of the droplet centers • fo realistic exper- 
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imental situations, the number of identical pion pairs in 
each single event is limited. We continue our investiga¬ 
tion here to find out appropriate measurable quantities 
that could be used to detect granular structure in HBT 
measurements. 

In our present investigation, we study the fluctuation 
between the single-event correlation function and the 
mixed-event correlation function and And that the distri¬ 
bution of the correlation function fluctuation / between 
the single-event and the mixed-event correlation func¬ 
tions can be a measurable quantities that could be used 
to probe the granular droplet structure. The width of the 
distribution is greater for a granular source than that for 
a single source, and the width increases as the droplet 
number decreases. The effect of the increase becomes 
less pronounced when the droplets have a collective ex¬ 
pansion. These changes of the widths can be quantified in 
terms of the root-mean-square fluctuation of /. The frms 
for a granular droplet source increases with the number 
of identical pion pairs in an event. The detection of the 
granular droplets becomes more favorable as the number 
of identical pairs increases. 

Our numerical simulated calculations indicate that the 
correlation function fluctuation leads to detectable differ¬ 
ences if the droplet number is small (less than or equal to 
16) and the number of identical pairs in each single event 
is of order 5 x 10"'’ or more. At RHIC energies, the multi¬ 
plicity of identical pion of event is about a few hundreds, 
and the number of pion-pairs in an event is about 10®. 
It will be of interest to see whether the correlation func¬ 
tion fluctuation can indeed be measured at RHIC. The 


situation become even more favorable for LHC collisions 
at higher energies where there can be a greater number 
of identical pion pairs. 

As with the development of many new experimental 
tools, progress is made by gradually increasing the com¬ 
plexity of one’s scopes of research and the areas of focus. 
In the present analysis, we have not considered the ab¬ 
sorption and the multiple scattering of the pions . 

Investigations on these effects of can be carried out in 
the future to see how they may modify the distribution 
of the fluctuations. We have also not considered the fluc¬ 
tuation of the overall size of emission source which clearly 
depends on the experimental selection and will require an 
investigation in conjunction with the experimental set-up 
and selections. An event multiplicity cut may be needed 
to reduce the fluctuation due to the sizes of the emitting 
source. Future investigations to refine the tool of HBT 
measurements for the detection of granular structure will 
be of great interest to probe the order of the quark-gluon 
plasma phase transition. 
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